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1  Abstract 


2  Four  surveys  of  airborne  expendable  bathythennograph  with  horizontal  spacing  of 

3  about  35  km  and  vertical  spacing  of  1  m  extending  from  the  surface  down  to  400  m  deep 

4  are  used  to  analyze  thermal  finestructurcs  and  their  seasonality  in  frontal  zones  of  the 

5  southern  Yellow  Sea  and  the  East  China  Sea.  Finestructure  characteristics  are  different  not 

6  only  among  fronts  but  also  along  the  same  front,  implying  different  mixing  mechanisms. 

7  Summer  thermocline  intrusions  with  thickness  from  few  to  40  meters  generated  by  the 

8  vertically-sheared  advection,  are  identified  along  the  southern  tongue  of  the  Cheju-Yangtze 

9  Front  (especially  south  of  Cheju  Island).  The  finestructures  south  of  the  Yangtze  Bank  (i.e., 

10  the  western  tip  of  the  southern  tongue)  produced  by  strong  along-frontal  currents  are  not  as 

11  rich  as  elsewhere  in  the  southern  tongue.  The  Cheju-Tsushima  Front  presents  mixed 

12  finestructures  due  to  confluent  currents  from  various  origins.  The  irregular-staircase 

13  finestructures  in  the  Kuroshio  region  (below  the  seasonal  thermocline)  driven  by  double- 

14  diffusive  mixing  show  seasonal  invariance  and  verticaFhorizontal  coherence.  The  strength 

15  of  mixing  related  to  finestructure  is  weaker  in  the  Kuroshio  region  than  in  the  Cheju- 

16  Tsushima  Front  or  south  of  Cheju  Island.  The  profiles  in  the  Tsushima  Warm  Current 

17  branching  area  show  large  (—50  m  thick)  irregular-staircase  structures  at  upper  230  m  depth, 

18  which  coincides  roughly  with  the  lower  boundary  of  maximum  salinity  layer.  The 

19  finestructure  at  deeper  230  m  depth  is  similar  to  that  in  the  Kuroshio  region.  The  possible 

20  mechanisms  for  generating  the  finestructures  are  also  discussed. 
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1.  Introduction 


2  Temperature  and  salinity  profiles  are  smooth  curves  (or  straight  lines)  or  ragged 

3  curves  (i.e.  a  variety  of  vertical  variations).  The  vertical  variation  with  scales  of  1-100  m  is 

4  termed  as  ‘finestructure’  and  that  with  scales  smaller  than  1  m  is  termed  as  ‘microstructure’ 

5  (Warren  and  Wunsch,  1981).  The  finestructure  is  observed  in  various  forms,  and  some  of 

6  them,  which  are  frequently  observed  and  detected,  are  named  as  follows:  (1)  staircase, 

7  also  called  steppy  or  stepwise  structure  (Schmitt  et  al.,  1987;  Zodiatis  and  Gasparini,  1996), 

8  when  temperature  and  salinity  profiles  show  alternative  sheets  [strong  temperature/salinity 

9  gradient  in  short  vertical  distance  such  as  di  to  d2  (~  I  m)  or  d3  to  d4  (~  2  m)  in  Fig.  la]  and 

10  layers  [almost  isothermal/isohaline  or  low  gradient  in  rather  longer  vertical  distance  such  as 

11  d2  to  d3  (~13  m)  in  Fig.  la];  (2)  thermocline/halocline  intrusion,  when  temperature  and 

12  salinity  profiles  display  multilayered  structures  at  the  thermocline/halocline  depths 

13  [Ruddick  and  Richards,  2003;  see  intrusive  structures  at  50-100  m  (50-90  m)  depth  in 

14  profile  2  (3)  in  Fig.  lb];  and  (3)  interleaving,  which  is  similar  to  the  thermocline/halocline 

15  intrusion  and  is  highlighted  by  large,  even  O  (100  km)  in  the  equatorial  Pacific,  horizontal 

16  continuity  of  multilayered  structures  (Richards  and  Banks,  2002;  Lee  and  Richards,  2004). 

17  In  addition,  there  are  mixed,  eroded  versions  of  these  forms,  and  any  forms  are  also 

18  possible  (Fig.  lc).  These  finestructures  are  common  anywhere  in  oceans  but  more  complex 

19  and  heterogeneous  around  fronts  and  current  confluence  zones  (Joyce,  1976;  Williams, 

20  1981;  Bianchi  et  al.,  1993). 
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Numerous  analyses  on  finestructures  have  been  conducted  in  the  tropics 

2  (Mcphaden,  1985;  Richards  and  Banks,  2002;  Lee  and  Richards,  2004),  the  Antarctic 

3  (Georgi,  1978,  1981;  Toole,  1981),  the  Atlantic  (Joyce,  1976;  Schmitt  and  Georgi,  1982; 

4  Schmitt  et  ah,  1986),  and  the  Pacific  (Gregg,  1977;  Peters  et  al.,  1991).  A  new  seismic 

5  reflection  methodology  produces  snapshots  of  finestructures  at  high  lateral  sampling  (<  10 

6  m)  in  the  Norwegian  Sea  (Nandi  et  al.,  2004).  These  analyses  enhance  the  knowledge  of 

7  mixing  processes  of  water  mass,  energy  transfonnation,  and  internal  wave  kinematics. 

8  However,  finestructures  or  mixing  processes  in  the  Yellow  Sea  and  the  East  China  Sea 

9  (YES)  have  been  less  studied,  although  the  YES  exhibits  various  thermohaline 

10  finestructures  induced  by  the  Kuroshio  intrusion,  enormous  river  discharge,  strong  tidal 

1 1  currents,  and  internal  waves  in  the  well-developed  shelves  along  with  trenches  (Chu  et  al., 

12  1997a,  b;  2005). 

13  Bao  et  al.  (1996)  analyzed  the  vertical  scales  and  spectral  characteristics  of  the 

14  finestructures  southwest  off  Kyushu  using  springtime  profiles  and  addressed  the 

15  finestructure  characteristics  with  respect  to  water  mass  distributions.  Lee  et  al.  (2003) 

16  reported  cross-frontal  thennal  intrusion  with  5-10  m  thickness  southeast  off  Cheju  Island 

17  and  several  high  or  low  salinity  cores  with  horizontal  scales  of  tens  kilometers  southwest  of 

18  the  Tsushima  Warm  Current  branching  region.  In  addition,  several  studies  showed 

19  thermohaline  finestructures  (Lie  et  al.,  1998;  Lie  et  al,  2000;  Lie  et  al.,  2003)  but  did  not 

20  provide  detail  information.  Park  and  Chu  (2006)  suggested  that  the  frontal  mixing  in  both 

21  temperature  and  salinity  fields  is  vigorous  in  the  region  connecting  to  the  Yangtze  Bank, 
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the  Cheju  Strait,  and  the  Korea  Strait.  Their  suggestion  is  based  on  the  analysis  of  the 

2  surface  and  subsurface  thermal  and  haline  fronts  throughout  YES  from  a  climatological 

3  dataset  [i.e.,  Generalized  Digital  Environmental  Model  (GDEM)]  in  tenns  of  frontal 

4  intensity,  water  mass  distribution,  and  temporal  evolutions  of  temperature  and  salinity 

5  across  the  fronts.  However,  finestructure  and  its  seasonality  were  not  investigated  due  to 

6  the  limitation  of  the  climatological  data. 

7  Recently,  new  thermal  features  were  discovered  using  four  airborne  expendable 

8  bathythermograph  (AXBT)  surveys  with  horizontal  spacing  of  35  km  and  vertical 

9  resolution  of  1  m  from  the  surface  to  400  m  depth  over  the  southern  Yellow  Sea  (YS)  and 

10  the  East  China  Sea  (ECS)  in  September  1992  and  February,  May,  and  September  1993 

11  (Furey  and  Bower,  2005;  Park  and  Chu,  2007a,  b).  Furey  and  Bower  (2005)  analyzed 

12  variability  of  mixed  layer  depth,  path  change  of  the  Kuroshio  associated  with  the  generation 

13  of  cold  eddies  on  canyons  northeast  off  Taiwan,  and  seasonal  evolution  of  thermal  fronts  in 

14  YES  and  the  Japan/East  Sea.  Park  and  Chu  (2007a)  identified  synoptic  distributions  of 

15  thermocline  and  thennal  surface  mixed  layer  [with  different  definition  compared  to  Furey 

16  and  Bower  (2005)],  as  well  as  seasonality  and  dominate  driving  mechanism  of  the  surface 

17  mixed  layer.  Park  and  Chu  (2007b)  further  investigated  the  synoptic  features  in/around 

18  thermal  fronts  and  the  cross-frontal  heat  fluxes  that  have  not  been  described  in  detail 

19  previously.  The  thermal  finestructures  in/around  the  fronts,  such  as  the  thennocline 

20  intrusions,  ragged  isotherms,  and  multilayered  temperature  inversions,  were  described  by 

2 1  Park  and  Chu  (2007b)  qualitatively  but  not  quantitatively. 
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Since  the  fronts  in  the  YES  are  linked  dynamically  to  the  current  system 

2  comprising  the  Kuroshio,  the  Taiwan  Warm  Current,  the  Tsushima  Warm  Current,  and  the 

3  Cheju  Warm  Current,  these  AXBT  data  with  sufficiently  large  horizontal  coverage  are 

4  useful  to  examine  the  finestructures  across  a  front  or  among  different  fronts  regarding 

5  surrounding  conditions  such  as  frontal  features,  water  masses,  and  the  current  system.  They 

6  are  also  invaluable  to  identify  the  seasonal  variability  of  the  finestructures  in  the  YES. 

7  The  goal  of  this  study  is  to  analyze  characteristics  of  thennal  finestructures, 

8  including  their  seasonality,  by  frontal  zones  in  the  YES  using  the  AXBT  data.  For  the  lack 

9  of  concurrent  salinity/velocity  data  and  time-series  data  for  analysis,  it  is  difficult  to 

10  elucidate  generation  mechanisms  of  the  finestructure  with  concrete  evidences.  Instead, 

1 1  using  the  AXBT  data,  CTD  data  (observed  at  different  time)  and  hydrographic  plots  in 

12  other  studies,  we  propose  plausible  generation  mechanisms  of  the  finestructures  and  related 

13  mixing  processes  in  connection  to  the  surrounding  thermal  features,  water  masses,  and 

14  circulation  systems. 

15  The  rest  of  the  paper  is  organized  as  follows.  Section  2  introduces  the  AXBT  data. 

16  Section  3  describes  finestructures  in  the  profiles  grouped  by  the  frontal  zones.  Section  4 

17  analyzes  ensemble-averaged  spectra  and  statistical  characteristics  of  groups  with  rich 

18  finestructures.  Section  5  examines  finestructures  of  each  profile  in  the  region  where  the 

19  Tsushima  Warm  Current  is  originated  from  the  Kuroshio.  Section  6  presents  conclusions. 
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1  2.  Data 


2  The  AXBT  data,  a  part  of  the  Master  Oceanographic  Observation  Data  Set 

3  (MOODS)  maintained  by  the  Naval  Oceanographic  Office  (NAVOCEANO),  Stennis  Space 

4  Center,  Mississippi,  were  obtained  under  the  approval  by  NAVOCEANO.  The  AXBT  data 

5  including  their  detail  information  have  not  been  open  to  the  public  yet.  Accordingly,  few 

6  references  have  been  available  although  the  AXBT  data  were  collected  in  the  early  1990s. 

7  In  1990s,  Sparton  AXBTs  were  widely  used  in  U.S.  Navy  research.  Their  thermal  time 

8  constant  0.1  s  or  less  and  their  vertical  resolution  is  15  cm  (Boyd  and  Linzell,  1993).  We 

9  obtained  the  edited  AXBT  data  with  1  meter  vertical  resolution.  The  editing  process 

10  includes  indentifying  recording  errors  and  outliers,  checking  duplicate  profdes,  checking 

1 1  depth  inversion  and  depth  duplication  for  individual  profdes,  checking  temperature  range, 

12  checking  large  temperature  inversions  and  gradients,  checking  standard  deviation, 

13  interpolation  to  standard  levels,  and  post  objective  analysis  checks  (Teague,  1986;  Jugan 

14  and  Beresford,  1992;  Boyer  and  Levitus,  1994;  Chu  et  al.,  1997c,  1998).  The  accuracy  of 

15  the  data  depends  on  the  accuracy  of  conversion  equations,  which  transfonn  frequency  into 

16  temperature  and  elapsed  fall  time  into  depth:  the  AXBT  does  not  measure  depth  directly. 

17  Customized  equations,  which  were  achieved  using  concurrent  CTD  data,  provide  the 

18  accuracy  of  temperature  data  of  0.1 3°C  and  that  of  depth  of  ±2%  of  depth  or  ±10  m, 

19  whichever  is  greater.  The  Naval  Research  Laboratory  (NRL)  Isis  System  determines  the 

20  AXBT  frequency  to  such  accuracy  that  the  resulting  temperature  accuracy  is  0.05°C  or 

21  better  (Boyd  and  Linzell,  1993). 
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1  The  AXBT  data  used  in  this  study  consist  of  1256  profiles  from  the  four  surveys 

2  (Fig.  2a):  18-29  September  1992  (named  as  9209),  4-14  February  1993  (as  9302),  5-14 

3  May  1993  (as  9305),  and  2-10  September  1993  (as  9309).  The  four  surveys  almost 

4  repeated  themselves  in  track  paths  with  slight  mismatches  and  were  completed  within  one 

5  week  except  for  some  profiles  near  the  Korea/Tsushima  Strait.  Note  that  the  seasonal 

6  variability  of  frontal  structures  can  be  identified.  For  convenience,  February,  May,  and 

7  September  represent  winter,  spring,  and  summer,  respectively. 

8  The  aspect  ratio  {H/L\  //:  vertical  scale;  L\  horizontal  scale)  of  individual  elements 

9  of  finestructures  is  10"4-10"3;  that  is,  a  finestructure  with  a  vertical  scale  of  few  meters  has 

10  approximately  a  horizontal  scale  of  a  few  to  tens  kilometers.  The  characteristic  time  scale  is 

1 1  clearly  associated  with  the  corresponding  spatial  scale  and  generally  few  weeks  (Fedorov, 

12  1978).  The  AXBT  data  with  the  horizontal  spacing  of  35  km  and  the  vertical  resolution  of  1 

13  m  are  therefore  capable  of  detecting  finestructures.  The  synoptic  features  identified  from 

14  these  data  are  realistic  under  the  assumption  that  the  oceanic  features  with  a  horizontal 

15  scale  larger  than  tens  kilometers  persist  for  at  least  a  week.  Here  the  vertical  scales  of  the 

16  finestructures  are  defined  as  2-100  m  because  ocean  features  with  vertical  scales  smaller 

17  than  1  m  cannot  be  identified  in  these  data.  We  do  not  examine  small-scale  fluctuations  in 

18  surface  and  bottom  layers,  although  they  are  related  to  turbulent  mixing.  Due  to  lack  of 

19  concurrent  salinity  and  velocity  data  and  time-series  data,  we  cannot  discuss  the  effect  of 

20  internal  waves  on  finestructures,  although  the  internal  waves  have  been  often  observed  in 
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1  the  YES  (Hsu  et  al.,  2000;  Han  et  al.,  2001).  The  thennal  (Incstructurcs  at  the  thermocline 

2  depth  in/around  the  thennal  frontal  zones  are  studied  here. 

3  3.  Two  Measures  for  Double  Diffusion 


4  There  are  two  measures  for  the  double-diffusive  mixing.  The  first  one  is  the  Turner 

5  angle  (Tu)  (Ruddick,  1983): 

6  Tu  = -tan  R  =  (add  /  dz)/(j3dS  /  dz),  (1) 


where  R  is  the  density  ratio,  and 


1  (  dp']  1  (  dp 

a  = - —  and  p  =  —  — 

p\86  )  P\dS 


9  are  the  coefficients  of  thennal  expansion  and  haline  contraction.  The  water  column  is 


10  gravitationally  stable  for  \Tu\<n!4,  gravitationally  unstable  for  \Tu\>  a  12 ,  double 


1 1  diffusive  for  n  /  4  <  Tu  <  n  /  2  ,  and  the  salt  fingering  for  -n  /  2  <  Tu  <  -k  /  4 . 

12  The  second  measure  is  the  Cox  number,  which  is  the  normalized  temperature 

13  gradient  variance  defined  by 


14 


CT 


Var(ATz) 
- \2  ’ 


(3) 


15  which  serves  as  a  non-dimensional  measure  of  finestructure  intensity  (Joyce,  1976).  Here, 

16  Var(ATz)  is  the  variance  of  vertical  temperature  gradient  in  a  segment  of  the  profile,  and 


17  ATZ  is  the  averaged  vertical  temperature  gradient  over  the  segment. 

18  (need  general ...) 
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1  4.  Detection  of  finestructures  in  frontal  regions 

2  The  temperature  profiles  are  locally  grouped  with  respect  to  their  vertical 

3  temperature  structures  related  to  the  thermal  fronts  (Fig.  2)  identified  from  these  data  (see 

4  Fig.  2b):  Cheju-Yangtze  Front  (CYF),  Cheju-Tsushima  Front  (CTF),  Tsushima  Front  (TF), 

5  and  Kuroshio  Front  (KF)  (Park  and  Chu,  2007b).  The  CYF  is  two-tongue-shaped,  i.e. 

6  slanted  S-shaped  (northern  and  southern  tongues).  The  CTF  occurs  along  the  southern  coast 

7  of  Korea  and  expands  between  the  Cheju  Strait  and  the  Korea/Tsushima  Strait.  The  TF  is  a 

8  branch  from  the  KF  (Park  and  Chu,  2006,  2007b). 

9  Seven  groups  are  also  identified  (Fig.  2a):  (1)  Kuroshio  Front  group  (KFG;  in  the 

10  Kuroshio  Front  but  the  Tsushima  Warm  Current  branching  region  is  not  included),  (2) 

11  Cheju-Tsushima  Front  group  (CTFG;  in  the  Cheju-Tsushima  Front),  (3)  Yellow  Sea 

12  Bottom  Cold  Water  group  (YCWG;  northwest  off  the  northern  tongue  of  the  CYF),  (4) 

13  west  of  Cheju  Island  group  (WCG;  in  the  northern  tongue  of  the  CYF),  (5)  south  of  Cheju 

14  Island  group  (SCG;  in  northern  part  of  the  southern  tongue  of  the  CYF),  (6)  Yangtze  Bank 

15  group  (YBG),  and  (7)  south  of  Yangtze  Bank  group  (SYBG).  The  last  two  groups  (YBG 

16  and  SYBG)  are  located  north  and  south  of  the  western  tip  of  the  southern  tongue  of  the 

17  CYF  with  consideration  of  seasonal  migration  of  the  southern  tongue  of  the  CYF.  The 

18  SYBG  also  shows  a  thennal  front  induced  by  the  Taiwan  Warm  Current  and/or  the  uplifted 

19  Kuroshio  northeast  of  Taiwan  merging  into  the  southern  tongue  of  the  CYF  in  summer 

20  (Park  and  Chu,  2007b).  The  forms  of  finestructures  are  diverse  among  these  groups. 
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a.KFG 


2  The  staircase  structure  (Fig.  la)  is  often  found  in  the  KFG  regardless  of  seasons 

3  and  depths  except  the  surface  mixed  layer  (Fig.  3)  with  irregular  layer  thickness  from  few 

4  to  40  meters  and  sheet  thickness  around  few  meters.  Some  staircase  structures  are  eroded 

5  by  mixing  process  and  seen  as  less-defined  sheets  and  layers.  These  staircase  structures  are 

6  called  irregular-staircase  (Gregg,  1975;  Ruddick  and  Gargett,  2003)  or  eroded-staircase 

7  structure.  This  irregular-staircase  structure  in  profiles  is  common  in  open  oceans. 

8  According  to  Bao  et  al.  (1996)’s  springtime  observation  at  the  Kuroshio  axis,  the  T-S 

9  relation  is  tight  in  irregularly-spaced  layers.  Profiles  of  the  vertical  temperature  (salinity) 

10  gradient  multiplied  by  the  thermal  expansion  (haline  contraction)  coefficient  (available  at 

11  http://www.nodc.noaa.gov/)  displays  a  good  coherence  between  temperature  and  salinity 

12  variations  (Fig.  4b). 

13  When  the  Turner  angle  is  calculated  from  the  CTD  data  with  1  m  vertical 

14  resolution,  this  is  just  a  gross  measure  of  the  stability,  not  an  actual  micro-scale  double 

15  diffusion  (Kennan  and  Lukas,  1996).  The  halocline,  from  a  salinity  maximum  (-150  m)  to 

16  a  salinity  minimum  (-500  m),  lies  in  the  salt  finger  regime  (i.e.  double  diffusion) 

17  particularly  in  a  weak  fingering  regime  (2<R<oo),  agreeing  with  the  irregular-staircase 

18  structure  (Fig.  4c).  In  the  KFG  60%  water  column  in  50-400  m  depth  lies  in  the  weak 

19  fingering  regime,  according  to  Fig.  5.  The  double-diffusive  diapycnal  mixing  may  cause 

20  this  irregular-staircase  structure. 
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1  The  Cox  number  was  calculated  for  each  profile  and  then  averaged  over  the  groups 

2  (Table  1;  the  criteria  of  dividing  the  upper  KFG  and  the  lower  KFG  will  be  explained  in 

3  section  4).  The  Cox  number  in  the  KFG  ranges  from  0.3  to  0.7,  which  is  an  agreeable 

4  measure  for  the  irregular-staircase  (Hayes  et  a/.,  1975;  Joyce,  1976). 

5  b.  CTFG 

6  The  CTFG  has  highly-complex/mixed  fincstructurcs  in  summer,  which  are  hardly 

7  specified  as  one  fonn  of  the  finestructure:  a  mixture  of  multiple  inversions,  thermocline 

8  intrusions  at  scales  of  few  to  30  meters  (9209  and  9309  in  Fig.  6).  These  finestructures  are 

9  overlapped  on  a  thermocline  but  do  not  change  considerably  a  main  trend  of  the 

10  thermocline.  They  are  richer  in  the  lower  part  than  in  the  upper  part  of  the  thennocline. 

1 1  Warm-noses  (inverted  structures  describing  maximum  temperature  observed  between  the 

12  surface  mixed  layer  and  the  top  of  a  strong  thermocline)  reported  by  Park  and  Chu  (2007b) 

13  are  seen  at  10-50  m  depth.  The  Cox  number  is  larger  than  1  in  9209  and  9309,  confirming 

14  higher  finestructure  intensity  than  that  in  the  KFG  (Table  1).  In  9209,  a  profile  displays 

15  ~9°C  near  the  bottom,  which  is  related  to  the  Korea  Strait  Bottom  Cold  Water,  the  water 

16  mass  passed  through  the  Korea/Tsushima  Strait  from  the  Japan/East  Sea  in  the  bottom  layer. 

17  In  9209  the  cold  water  extends  west  of  129. 5°E  along  the  trench  of  the  Korea/Tsushima 

18  Strait,  while  in  9309  does  not  (see  Fig.  1 1  in  Park  and  Chu,  2007b). 

19  c.  YCWG  and  WCG 

20  The  YCWG  shows  smooth  curves  or  straight  lines  (i.e.  lack  of  finestructures)  but  a 

21  strong  seasonal  variability  except  in  the  bottom  layer  (Fig.  7).  A  colder  layer  formed  above 
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1  a  homogenous  bottom  cold  layer  (i.e.  the  Yellow  Sea  Bottom  Cold  Water)  exists  in  9305. 

2  No  prominent  flnestructure  except  this  colder  layer  is  detected  in  warm  seasons  since  the 

3  Cox  numbe  is  quite  low  (~0.3)  (Table  1). 

4  The  WCG  is  similar  to  the  YCWG  (Fig.  8)  also  with  low  Cox  number  (0.2-0. 3). 

5  However,  the  WCG-9302  was  separated  into  two  groups  with  the  occurrence  of  the 

6  northern  tongue  of  the  CYF.  Vertical  temperature-salinity-density  distributions  across  the 

7  northern  tongue  (clear  coherent  haline  and  thermal  fronts)  show  the  existence  of  a  density- 

8  compensated  front  (Lee  et  al.,  2003;  Son  et  al.,  2003).  Wintertime  wind  stirring  suppresses 

9  the  flnestructure  development  around  the  northern  tongue  despite  the  homogenous  density 

10  field.  Nonetheless,  the  inversion  layer  at  40-70  m  depth  with  small  temperature 

1 1  fluctuations  in  WCG-9302  implies  stronger  cross-frontal  mixing  and/or  advection  than 

12  wind  stirring  (Park  and  Chu,  2007a).  The  colder  layer  fonned  above  the  homogenous 

13  bottom  cold  layer  in  warm  seasons  is  not  as  rich  as  that  in  the  YCWG. 

14  d.  SCG 

15  The  SCG  is  characterized  by  intensive  summer  thennocline  intrusions,  owing  to 

16  the  mixing  between  the  modified  Yellow  Sea  Bottom  Cold  Water,  i.e.,  the  water  mass  of 

17  10-15°C  over  the  Yangtze  Bank  extended  southward  from  the  Yellow  Sea  Bottom  Cold 

18  Water  and  mixed  with  cold  Chinese  coastal  water  (fsobe,  1999;  Furey  and  Bower,  2005; 

19  Park  and  Chu,  2006,  2007b),  and  the  Tsushima  Wann  Current  water  (9209  and  9309  in  Fig. 

20  9).  The  intrusions  are  developed  at  20-100  m  depth,  beneath  the  strong  upper  seasonal 

21  thermoclines,  with  horizontal  scale  of  tens  kilometers  and  vertical  scale  of  a  few  to  40 
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1  meters  (Park  and  Chu,  2007b).  Especially,  in  the  SCG-9309  temperature  fluctuates  with 

2  amplitude  of  0.2-2°C  in  a  range  of  17-20°C  at  25-80  m  depth,  which  implies  mixing 

3  between  the  strong  upper  seasonal  thermocline  and  the  bottom  layer.  The  fluctuation  in  this 

4  intermediate  layer  seems  to  be  linked  to  the  homogenous  water  of  ~  1 8°C  that  splits  locally 

5  the  combined  front  of  the  CYF  and  the  TF  in  9309  (Park  and  Chu,  2007b).  The  intensive 

6  summer  thermocline  intrusions  lead  to  the  largest  Cox  number  (4-8)  within  all  groups 

7  (Table  1). 

8  According  to  the  temperature  and  salinity  observations  in  June  1999  (Fee  et  al., 

9  2003),  halocline  intrusions  are  multilayered  at  rather  smaller  vertical  scale  (—10  m)  in  the 

10  western  part  of  the  SCG  (125. 5-126. 5°E,  ~32.5°N)  in  comparison  with  thennocline 

1 1  intrusions.  They  interpreted  that  these  intrusions  may  be  related  to  double-diffusive  mixing 

12  occurring  along  the  isopycnal  surface.  In  their  temperature  and  salinity  observations  in 

13  September  1998,  haline  features  resemble  thermal  features  beneath  the  strong  upper 

14  thermoclines/haloclines,  and  salinity  varies  in-phase  with  temperature  but  rather  complex. 

15  On  the  base  of  these  thermohaline  features,  the  intrusion  is  generated  also  by  the  vertically- 

16  sheared  advection;  that  is,  the  Yangtze  Diluted  Water  and  the  cold  shelf  water  on  the 

17  Yangtze  Bank  are  advected  seaward  in  the  intermediate  layer,  while  the  Tsushima  Warm 

18  Current  water  is  advected  shoreward  in  the  bottom  layer  (Fee  et  al.,  2003;  Park  and  Chu, 

19  2007b).  The  thermocline  intrusions  are  also  detected  in  9305  but  not  to  greater  depths. 

20  Temperature  in  the  bottom  mixed  layer  is  the  coldest  in  9305,  indicating  southeastward 
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1  migration  of  the  southern  tongue  of  the  CYF  (Furey  and  Bower,  2005;  Park  and  Chu, 

2  2007b,  2006). 

3  e.  YBG  and  SYBG 

4  The  YBG  shows  evident  thermal  seasonality  in  the  bottom  mixed  layer  because  of 

5  shallow  depth  and  strong  tidal  mixing  (Park  and  Chu,  2007a):  9-12°C  in  9302,  10-12. 5°C 

6  in  9305,  and  22-25°C  (<23°C)  in  9209  (9309)  (Fig.  10).  Two  discernable  groups  of  the 

7  profiles  indicate  the  southern  tongue  of  the  CYF  except  in  9302.  No  prominent 

8  finestructure  is  detected.  In  addition,  the  Cox  number  is  the  lowest  in  the  groups. 

9  The  SYBG  also  reveals  strong  horizontal  thennal  gradient  (Fig.  1 1)  particularly  in 

10  summer,  which  is  driven  by  the  Taiwan  Warm  Current  and/or  Kuroshio 

1 1  meandering/frontal  eddy-like  features  (Furey  and  Bower,  2005;  Park  and  Chu,  2007b).  The 

12  finestructure  in  the  SYBG,  however,  is  not  as  evident  as  that  in  the  SCG,  implying  weak 

13  cross-frontal  mixing.  When  a  strong  current  (i.e.,  the  Taiwan  Warm  Current  and/or 

14  Kuroshio  meander)  plays  a  major  role  in  frontal  genesis,  a  cross-frontal  (along-frontal) 

15  current  would  be  weak  (strong)  and  consequently  frontal  mixing  would  be  weak  as  in  the 

16  SYBG.  On  the  other  hand,  when  a  weak  current  develops  along  the  front,  a  cross-frontal 

17  current  would  not  be  negligible  and  consequently  may  induce  cross-frontal  mixing  by 

18  thermohaline  intrusions  as  in  the  SCG  (Tomczak  and  Godfrey,  2003).  Particularly,  in  9209 

19  no  finestructure  in  the  thermocline  is  attributable  to  strong  winds  at  that  time.  Thick  surface 

20  mixed  layer  pushes  down  the  thermocline,  and  inhibits  development  of  finestructures  in  the 

2 1  thermocline  (Park  and  Chu,  2007a). 
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1  5.  Spectral  and  statistical  characteristics 

2  Characteristics  of  the  finestructures  detected  in/around  the  thermal  fronts  can  be 

3  represented  by  vertical  temperature  gradient  or  temperature  fluctuation  (i.e.  a  high-pass- 

4  filtered  or  mean-removed  temperature):  for  instance,  intervals  between  oscillations, 

5  amplitudes  of  oscillations,  and  signs  of  oscillations  in  the  vertical  temperature  gradient 

6  profile  reveal  the  size  and  the  form  of  finestructure  (see  correspondence  of  the  temperature 

7  profile  to  its  vertical  temperature  gradient  in  Fig.  la).  However,  for  data  collected  over  the 

8  shelf  shallower  than  100  m,  it  is  inappropriate  to  apply  a  high-pass-filter  to  such  short 

9  profiles  for  removing  larger  (than  the  scale  of  finestructure)  fluctuations,  and  therefore  the 

10  vertical  temperature  fluctuation  is  not  available  for  these  data.  In  addition,  without  time 

1 1  series  of  the  profiles  we  cannot  calculate  the  mean-removed  temperature  profile 

12  accordingly.  Thus,  the  vertical  temperature  gradient  calculated  at  every  vertical  interval  (1 

13  m)  was  used  to  analyze  the  finestructure  characteristics. 

14  We  conducted  spectral  and  statistical  analyses  of  the  vertical  temperature 

15  gradient  ATZ  of  locally-grouped  profiles  with  rich  finestructures  (surface  and  bottom  mixed 

16  layers  were  excluded  in  calculations):  the  CTFG,  the  SCG,  and  the  KFG.  Since  the  KFG 

17  contains  the  surface  mixed  layer  and  the  seasonal  thermocline  at  depths  shallower  than 

18  -150  m,  its  profiles  were  divided  into  two  segments:  upper  (upper-KFG,  shallower  than 

19  150  m)  and  lower  (lower-KFG,  151-400  m).  Winter  profiles  with  poor  finestructures  were 

20  excluded  in  the  calculations  (except  lower-KFG). 
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1  First,  A Tz  spectra  of  typical  profiles  of  the  groups  were  calculated  to  estimate  the 

2  dominant  size  of  finestructure  (Fig.  12).  Second,  ensemble-averaged  probability 

3  distribution  functions  (PDFs)  of  A Tz  were  calculated  (Fig.  13).  Third,  a  diagram  of 

4  skewness  vs.  interquartile-range  was  plotted  (Fig.  14a).  For  the  temperature  gradient,  the 

5  positive  skewness  is  related  to  the  staircase  structure;  the  more  staircase-like  a  profile  is, 

6  the  more  probable  low  gradient  values  (i.e.  layer)  will  be  than  high  gradient  values  (i.e. 

7  sheet)  because  the  low  gradient  occupies  a  larger  portion  of  the  profile  (McPhaden,  1985). 

8  The  interquartile-range  is  the  difference  between  the  75th  and  the  25th  percentiles  of  data. 

9  The  large  interquartile-range,  i.e.  the  gradient  spreads  widely,  indicates  the  development  of 

10  any  forms  of  finestructure  except  the  staircase  structure  (the  interquartile-range  of  the 

1 1  profile  with  rich  staircase  structures  is  small,  because  the  low  gradient  (i.e.  layer)  occupies 

12  a  larger  portion  of  the  profile).  Most  summer  profiles  present  large  interquartile-range  (Fig. 

13  14a).  Fig.  14a  shows  four  zones  (zone-1  to  zone-4)  according  to  finestructure 

14  characteristics.  The  zoning  histogram  shows  the  profile  occurrence  percentage  of  each 

15  zone  (Fig.  14b).  More  explanations  of  the  zones  will  be  given  later.  Fourth,  a  diagram  of 

16  the  Cox  number  vs.  gradient  ratio  was  plotted  (Fig.  15a).  The  gradient  ratio  is  defined  as 

17  the  occurrence  of  A Tz  <0(i.e.  inverted  or  isothennal  layer)  divided  by  that  of  A Tz  >0 

18  (i.e.  layer  where  temperature  decreases  with  depth)  over  the  segment.  For  example,  when 

19  one  third  of  all  A Tz  calculated  at  every  vertical  interval  over  the  segment  is  non-positive 

20  and  the  rest  is  positive,  the  gradient  ratio  is  0.5.  The  larger  the  Cox  number  (or  the  gradient 
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1  ratio),  the  higher  the  potential  for  mixing  is.  The  intensity  and  the  fonn  of  fincstructurc  can 

2  be  represented  by  the  Cox  number-gradient  ratio  diagram  by  three  zones,  zone-  a  to  zone- 

3  y ,  (Fig.  15a)  and  corresponding  occurrence  (Fig.  15b). 

4  According  to  our  calculation  for  the  four  groups,  the  zone-  a  exhibits  the 

5  irregular-staircase  with  the  Cox  number  <0.7  and  the  gradient  ratio  <~0.1;  two  critical 

6  values  (Cox  number  =  0.7;  gradient  ratio  =  0.1)  are  indicated  in  Fig.  15a.  Our  range  of  the 

7  Cox  number  for  the  irregular-staircase  is  consistent  to  several  earlier  studies  (Hayes  et  ah, 

8  1975;  Joyce,  1976).  The  profiles  of  zone- a  generally  have  low  interquartile-range  (<0.1) 

9  and  high  skewness  (~2);  two  critical  values  (interquartile-range  =  0.1;  skewness  =  2)  are 

10  indicated  in  Fig.  14a.  Additionally,  some  of  these  profiles  show  even  lower  interquartile- 

1 1  range,  so  that  a  lower  critical  value  (interquartile-range  =  0.05,  marked  in  Fig.  14a)  divides 

12  the  zone-1  and  the  zone-2.  More  staircase-like  profiles,  i.e.  thicker  layer  of  the  stair,  would 

13  fall  in  zone  3.  The  zone-  ft  shows  the  thermocline  intrusion  with  the  Cox  number  >  1  and 

14  the  gradient  ratio  ~1.  The  profiles  of  zone-  fi  generally  have  higher  interquartile-range  and 

15  lower  skewness  than  the  profiles  of  zone-<z  (irregular-staircase)  (Fig.  14a;  zone -4).  The 

16  zone-y  is  an  intermediate  zone  between  zone- a  and  zone- /?  (with  partially  overlapping 

17  with  them).  Profiles  have  less  typical  characteristics  in  the  zone-  y  than  in  either  zone- a 

1 8  or  zone-  /? . 

19  A  well-defined  staircase  structure  is  expected  to  have  the  Cox  number  >1  and  the 

20  gradient  ratio  »  1 .  This  was  not  detected  in  the  AXBT  data,  though.  Even  in  zone-4,  a 
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1  profile  could  have  the  zone- <2  features  if  the  mean  gradient  of  thennocline  is  large  (e.g. 

2  some  profiles  in  CTFG).  The  two  diagrams  should  therefore  be  considered  together. 

3  a.  Upper  Segment  of  KFG 

4  The  spectra  of  the  vertical  temperature  gradient  in  the  upper  segment  show  a 

5  distinct  seasonal  variability  in  all  vertical  wavenumbers  (Fig.  12a).  In  the  summer  a 

6  dominant  size  of  the  finestructure  is  10-20  m  with  the  peak  at  0.05-0.1  cycles  per  meter 

7  (cpm)  (Fig.  12a).  For  A Tz  >  0.05°C  m  ,  the  probability  of  the  vertical  temperature  gradient 

8  enhances  from  winter  to  summer  (Fig.  13a).  Since  this  high  gradient  is  consistent  with  the 

9  gradient  of  the  sheets  (~0. 1°C  m'1)  in  the  KFG,  the  development  of  the  staircase  structure  is 

10  responsible  for  the  gradient  enhancement.  In  the  summer,  the  zone -2  occurs  more  often  in 

1 1  the  upper-KFG  than  in  the  lower-KFG  (Fig.  14b),  which  causes  larger  interquartile-range  in 

12  the  upper-KFG  than  in  the  lower-KFG.  This  implies  that  the  finestructure  scale  and/or  form 

13  tend  to  have  more  variability  in  the  upper-KFG  than  in  the  lower-KFG.  For  the  zone-3, 

14  profiles  are  more  staircase-like  (thicker  layer  in  the  stair)  in  9309  than  in  9209  (Fig.  14b). 

15  In  the  spring  there  is  no  significant  difference  in  the  finestructure  between  the  upper  and  the 

16  lower  segments  in  terms  of  spectra  and  statistics,  owing  to  the  absence  of  seasonal 

17  thennocline  (9305  in  Figs.  12a,  13a;  9305  in  Figs.  12b,  13b;  9305  in  Figs.  14b,  15b). 

1 8  b.  Lower  Segment  of  KFG 

19  In  contrast  to  the  upper-KFG,  the  spectra  of  the  vertical  temperature  gradient  in  the 

20  lower  segment  are  almost  seasonally  invariant  regardless  of  seasonal  temperature 

21  variability  (Figs.  3,  12b).  The  seasonal  invariance  of  the  finestructure  is  also  confirmed  by 
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1  the  PDFs  (see  four  very  similar  curves  in  Fig.  13b)  and  the  two  zoning  histograms  (see 

2  resemblance  of  bars  to  each  survey  in  Figs.  14b,  15b).  In  the  winter,  the  surface  mixed  layer 

3  deepening  affects  the  finestructure  slightly  through  increasing  zone-3  (zone- y )  in  9202 

4  compared  with  the  other  surveys  as  shown  in  Fig.  14b  (Fig.  15b). 

5  A  significant  occurrence  percentage,  60—80%,  of  zone-1  (zone-a)  in  Fig.  14b 

6  (Fig.  15b)  demonstrates  the  dominance  of  staircase  structure.  The  slope  of  the  spectra  is 

7  gentle  (~10"“)  in  high  vertical  wavenumbers  (>10‘  cpm)  (Fig.  12b).  The  Cox  number  of  the 

8  lower-KFG  is  lower  than  that  of  the  CTFG  and  the  SCG  (Fig.  15a).  These  spectral  and 

9  statistical  features  suggest  that  mixing  in  the  lower-KFG  is  more  frequent  at  vertical  scales 

10  <10  m  than  at  vertical  scales  >10  m,  but  is  not  as  vigorous  as  in  the  CTFG  and  the  SCG. 

1 1  c.  CTFG 

12  The  energy  level  of  the  spectra  of  the  vertical  temperature  gradient  is  1 0 1  °- 1 0 1 5 

13  times  higher  in  the  summer  than  in  the  spring  (Fig.  12c).  In  the  summer,  the  PDFs  of 

14  vertical  temperature  gradient  display  a  low  peak  with  a  long  irregular  tail  toward  the  large 

15  vertical  gradient  (Fig.  13c),  evidencing  the  complex/mixed  finestructures  at  scales  of  a  few 

16  to  30  meters  mentioned  in  section  3b.  This  PDF  is  contrasted  to  a  high  peak  with  a  short 

17  smooth  tail  of  all-season  lower  segment  of  the  KFG  (Fig.  13b).  A  dominant  mixing 

18  mechanism  could  therefore  be  different  between  the  CTF  and  the  KF.  Note  that  the  CTF  is 

19  lack  of  one-zone-dominancy  unlike  the  other  groups  (Figs.  14b,  15b).  The  region  around 

20  the  CTF  is  a  confluence  area  of  various  water  masses:  the  water  mass  transported  by  the 

21  Tsushima  Warm  Current,  the  water  mass  passed  through  the  Cheju  Strait,  the  water  mass 
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1  along  the  south  coast  of  Korea,  and  the  water  mass  passed  through  the  Korea/Tsushima 

2  Strait  in  the  bottom  layer.  The  mixing  among  those  water  masses  is  complicated  and 

3  mixing  scale  varies  widely. 

4  d.  SCG 

5  The  PDFs  of  vertical  temperature  gradient  in  the  SCG  are  similar  as  in  the  CTFG, 

6  i.e.  the  low  peak  with  the  long  irregular  tail  toward  to  the  large  gradient  (Fig.  13d). 

7  However,  the  PDFs  in  the  SCG  show  higher  percentage  in  the  negative  to  zero  gradient 

8  range  than  those  in  the  CTFG;  especially  in  the  summer  the  probability  of  the  vertical 

9  gradient  in  this  range  increases  by  -45%.  Similarly  the  gradient  ratio  is  nearly  one,  highest 

10  in  the  YES  (Fig.  15a).  The  Cox  number  is  also  the  highest.  Accordingly,  zone-/?  occurs 

11  more  often  in  the  SCG  than  in  the  CTFG  (Fig.  15b).  More  than  60%  of  the  SCG  profiles 

12  belong  to  zone -4,  unlike  the  other  groups  (Fig.  14b).  These  features  are  associated  with 

13  thicker  uniform  and/or  inversion  layers  in  the  SCG,  which  result  from  the  temperature 

14  fluctuation  between  the  strong  upper  seasonal  thermoclines  and  the  bottom  layer  as 

15  mentioned  in  section  3d  and  yield  more  vigorous  mixing  than  the  CTFG.  Scale  and 

16  magnitude  of  the  vertical  temperature  fluctuation  are  also  greater  in  the  SCG  than  in  the 

17  CTFG  (Figs.  6,  9).  In  the  SCG  the  peaks  of  spectra  are  located  at  vertical  wavenumbers 

18  <5xlCT2cpm,  i.e.  vertical  scales  >  20  m  (Fig.  12d).  These  finestructure  characteristics  of 

19  the  SCG  are  consistent  with  the  greatest  cross-frontal  heat  flux  and  lateral  eddy  diffusivity 

20  in  the  YES  (Park  and  Chu,  2007b). 
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1  6.  Finestructures  in  the  Tsushima  Warm  Current  branching  region 

2  a.  Horizontal  incoherency 

3  The  Tsushima  Warm  Current  branching  occurs  generally  in  127-128°E  and  30°N. 

4  The  warm  saline  water  detached  from  the  Kuroshio  flows  onto  the  ECS  shelf,  despite 

5  seasonal  variations  (Huh,  1982;  Lie  et  al.,  1998;  Hsueh,  2000;  Furey  and  Bower,  2005)  (Fig. 

6  16).  In  the  wann  seasons  the  Tsushima  Warm  Current  connects  to  the  southern  tongue  of 

7  the  CYF  and  is  contributed  by  the  Taiwan  Wann  Current  (Isobe,  1999;  Ichikawa  and 

8  Beardsley,  2002;  Furey  and  Bower,  2005;  Park  and  Chu,  2007b).  Just  apart  from  the 

9  branching  location,  meanders  and  eddies  often  occur  (Lie  et  al.,  1998;  Nakamura  et  al., 

10  2003)  and  bottom  relief  is  complex  (Fig.  16).  Mesoscale  eddy  stirring  can  create 

1 1  finestructures,  which  are  characterized  by  lack  of  horizontal  coherency  (Ferrari  and  Polzin, 

12  2005).  For  these  reasons  the  temperature  profiles  in  the  Tsushima  Warm  Current  branching 

13  region  are  horizontally  less  coherent  than  those  in  the  KFG  (Fig.  17;  Fig.  2  by  Yang  et  al., 

14  2004).  Here,  for  convenience  the  branching  region  is  assumed  wider  to  encompass  the 

15  branching  location.  In  particular,  this  horizontal  incoherency  is  obvious  in  9302  and  9305. 

16  Hence,  finestructure  characteristics  here  need  to  be  analyzed  by  individual  profiles,  rather 

17  than  by  their  ensemble-averages. 

1 8  b.  Comparison  of  profiles  on  and  off  the  branching  location 

19  We  analyzed  individual  profiles  in  three  stations  (Figs.  18,  19).  One  of  them  was 

20  approximately  in  the  branching  location  (Stn-A;  128. 0°E,  29.7°N),  and  two  others  were  a 

21  little  north  (Stn-B;  127. 6°E,  30.8°N)  and  south  (Stn-C;  128. 6°E,  29.0°N)  of  the  branching 
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1  location  in  order  to  compare  fincstruc lures  on  and  off  the  branching  location  (Fig.  16).  The 

2  south  (north)  station  of  the  branching  location  presents  a  profile  before  (after)  the  branching. 

3 

4  At  Stn-B,  the  subsurface  Kuroshio  water  (<21°C,  >24.4  <j0 ;  Oka  and  Kawabe,  1998) 

5  intrudes  onto  the  shelf  at  depths  deeper  than  100  m  all  the  year  and  sustains  temperature  of 

6  16-17°C  in  the  bottom  layer  (Stn-B  in  Fig.  18a).  For  9305  and  9309  there  is  a  big  jump  in 

7  the  intermediate  layer  (-75  m  depth)  in  the  vertical  temperature  gradient  and  its  cumulative 

8  variance  (cumulate  variance  from  the  surface  to  every  depth  of  the  profile)  (Stn-B  in  Figs. 

9  18b,  c).  This  jump  implies  that  different  water  masses  meet  and  mix.  At  Stn-C,  on  the  other 

10  hand,  weak  gradient  exhibits  (Stn-C  in  Fig.  18b).  The  cumulative  variance  of  the  gradient  at 

1 1  Stn-C  is  less  than  one  sixth  of  that  at  Stn-B,  as  compared  at  130  m  depth  (Stn-C  in  Fig.  18c). 

12  The  energy  spectra  of  the  vertical  temperature  gradient  at  Stn-C  (except  9302)  are  also  101- 

13  10"  times  lower  than  those  at  Stn-B,  especially  much  lower  in  the  low  (<10"  cpm)  vertical 

14  wavenumbers  (Fig.  19).  These  spectral  and  statistical  characteristics  imply  poor 

15  finestructures  at  a  vertical  scales  >10  m  at  Stn-C  in  comparison  to  that  at  Stn-B. 

16  In  the  summer,  large  (10-50  m  thick)  irregular-staircase  structures  are  strong  with 

17  big  jumps  in  the  gradient  and  the  large  cumulative  variance  of  the  gradient  in  the  upper  part 

18  (shallower  than  230  m)  of  Stn-A  and  Stn-B  (Figs.  18a-c)  but  not  of  Stn-C.  The  PDFs  of 

19  vertical  temperature  gradient  of  Stn-A  and  Stn-B  have  the  second  peak  at  0.25°C  m"1 

20  (corresponding  to  the  sheets)  (Fig.  18d).  In  contrast,  the  large  irregular-staircase  structure  is 
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1  comparatively  weak  in  the  deeper  part  (deeper  than  230  m)  of  Stn-A,  which  is  similar  to 

2  Stn-C. 

3  The  energy  spectra  of  the  vertical  gradient  resemble  in  shape  and  slope  between 

4  Stn-A  and  Stn-C  in  the  high  vertical  wavenumbers  ( >  1 0" 1  cpm),  and  are  higher  in  Stn-A 

5  than  in  Stn-C  in  the  low  vertical  wavenumbers  (<  1  O'1  cpm)  (9209  and  9309  in  Fig.  19, 

6  although  9309  does  not  show  these  feature  as  clearly  as  9209).  Peaks  at  -0.02  cpm  (0.05- 

7  0. 1  cpm)  seem  to  be  related  to  the  size  of  layers  (sheets)  (9209  and  9309  in  Fig.  19-Stn  A). 

8  c.  Generation  mechanism 

9  Since  the  cold  Kuroshio  subsurface  water  is  uplifted  at  Stn-A,  temperature  of  Stn- 

10  A  is  colder  than  that  of  Stn-C  at  depths  deeper  than  100  m  (Fig.  18a).  This  uplifted  cold 

1 1  water  and  the  cold  shelf  water  cause  strong  horizontal  temperature  gradient  against  the 

12  warm  Kuroshio  surface  water  (>21°C)  (Oka  and  Kawabe,  1998;  Su  and  Wang,  1994),  and 

13  this  gradient  accordingly  induces  the  large  irregular-staircase  structure  in  the  upper  part  of 

14  Stn-A.  The  development  of  the  large  irregular-staircase  structure  is  confirmed  by  the 

15  difference  of  Stn-A  and  Stn-C  in  the  spectra  and  statistics  mentioned  in  Section  6b. 

16  This  strong  horizontal  gradient  seems  to  be  formed  also  in  salinity  field  and  then 

17  develop  an  irregular  halocline  structure,  although  we  do  not  have  concurrent  salinity  data. 

18  According  to  the  springtime  temperature  and  salinity  profiles  in  Bao  et  al.  (1996,  see  their 

19  Fig.  2),  which  are  close  to  Stn-A,  salinity  varies  generally  in  phase  with  temperature  and 

20  presents  rich  multilayered  structures  in  the  profiles  at  shallower  than  -200  m  depth.  Salinity 

21  (temperature)  varies  in  the  range  of  0.5-1  psu  (1-2°C)  in  a  sheet  with  thickness  <10  m.  In 
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addition,  the  vertical  salinity  distributions  in  the  branching  region  observed  in  May  1995 

2  present  a  variety  of  lincstructures,  such  as  low  or  high  cores  and  multilayered  structures,  in 

3  10-50  km  long  and  a  few  to  tens  meters  thick  according  to  Lie  et  al.  (1998,  see  their  Fig.  6). 

4  It  is  noted  that  the  lower  limit  of  these  finestructurcs  (-200  m  depth)  generally  coincides 

5  with  the  depth  of  maximum  salinity  zone  (-34.7  psu)  and  seasonal  pycnocline,  based  on 

6  their  observations  and  GDEM  data.  Fig.  4  also  presents  this  feature,  although  the  profile  is 

7  in  the  KFG.  No  prominent  finestructure  except  the  staircase  is  below  200  m,  like  in  Stn-C 

8  and  the  lower  segment  of  the  KFG  (Fig.  18c):  the  finestructure  develops  down  to  minimum 

9  salinity  zone  (-700  m  depth),  and  it  is  nearly  seen  beyond  the  zone  (Fig.  4). 

10  The  thermal  and  density  fronts  generated  by  the  Tsushima  Warm  Current 

1 1  branching  are  coherent  in  vertical  distributions  of  temperature  and  density,  whereas  the 

12  corresponding  haline  front,  narrower  than  the  thermal  or  density  front,  is  located  close  to 

13  the  shoreward  boundary  of  the  thermal  or  density  front.  Such  haline  finestructures  as 

14  mentioned  above  are  enriched  in  the  seaward  side  away  from  the  haline  front,  but  still 

15  inside  the  thermal/density  frontal  zone,  following  isopycnal  surfaces  (Lee  et  al.,  2003). 

16  Taking  these  features  into  account,  the  thermohaline  structures  tend  to  be  compensated  in 

17  the  branching  region,  consistent  to  Bao  et  al.,  (1996).  The  effect  of  temperature  on  density, 

18  however,  seems  to  be  greater  than  that  of  salinity,  although  richer  finestructures  in  their 

19  studies  are  detected  in  the  salinity  field  and  salinity  compensates  or  dominates  density 

20  locally  and  often  above  the  maximum  salinity  zone. 
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7.  Conclusions 


2  The  four  AXBT  surveys  were  used  to  analyze  thennal  finestructures  and  their 

3  seasonality  in  frontal  zones  of  the  southern  Yellow  Sea  and  the  East  China  Sea.  The 

4  possible  mechanisms  for  generating  these  finestructures  were  discussed.  For  the  lack  of 

5  concurrent  salinity/velocity  and  time  series  data,  only  limited  analyses  on  temperature  were 

6  possible.  The  four  surveys  conducted  in  one  year  are  not  quantitatively  sufficient  to 

7  describe  the  finestructure  seasonality  with  robust  statistics.  However,  with  a  good  synoptic 

8  view  of  the  AXBT  data,  we  demonstrated  that  the  finestructure  characteristics  are  different 

9  not  only  among  the  fronts  but  also  along  the  front,  suggesting  that  different  mixing 

10  mechanisms  are  related.  The  results  are  summarized  as  follows: 

11  (1)  The  thennal  finestructure  in  the  Kuroshio  Front  is  characterized  by  the 

12  irregular-staircase  consisting  of  alternative  sheets  (a  few  meters  thick  with  a  gradient  of 

13  ~0.1°C  m'1)  and  layers  (a  few  to  40  m  thick).  The  lower  (151-400  m)  segment  of  the 

14  Kuroshio  Front  Group  reveals  seasonal  invariance  in  finestructure,  unlike  the  upper 

15  segment  (shallower  than  150  m).  The  winter  surface  mixed  layer  deepening,  however, 

16  makes  the  layer  of  the  irregular-staircase  rather  thicker.  The  vertical  temperature  gradient  in 

17  the  range  of  0<  A Tz  <0.05  °C  m'1  accounts  for  70%  of  its  distribution  probability  in  all 

18  seasons,  implying  that  single  mixing  mechanism  (probably  double  diffusive  diapycnal 

19  mixing)  is  dominant  in  generating  the  finestructure.  The  Cox  number  is  lower  in  the 

20  Kuroshio  Front  than  in  the  Cheju-Tsushima  Front  or  the  northeastern  part  of  the  southern 

21  tongue  of  the  Cheju-Yangtze  Front,  and  the  scale  of  finestructure  is  comparatively  smaller. 
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1  (2)  The  temperature  profiles  in/around  the  Cheju-Tsushima  Front  reveal  highly 

2  complex/mixed  finestructures  at  the  scales  of  a  few  to  30  meters  in  the  lower  part  of  the 

3  thermocline  in  the  summer.  This  finestructure  characteristic  is  confirmed  by  the  low  peak 

4  with  a  long  irregular  tail  in  the  PDF  of  the  vertical  temperature  gradient  as  well  as  the  lack 

5  of  one-zone-dominancy  in  the  statistical  classifications  by  finestructure  characteristics. 

6  Since  the  Cheju-Tsushima  Front  region  is  the  confluence  area  in  the  YES  crowded  with 

7  water  masses  from  different  origins,  mixing  among  these  water  masses  is  likely  to  occur  in 

8  many  ways  and  to  yield  such  complex/mixed  finestructures. 

9  (3)  The  northern  tongue  of  the  Cheju-Yangtze  Front  shows  no  prominent 

10  finestructures  except  the  colder  layer  formed  above  the  homogenous  bottom  cold  layer  in 

1 1  warm  seasons.  Occurrence  of  winter  inversion  layer  at  40-70  m  depth  with  small 

12  temperature  fluctuations  indicates  sufficiently  strong  cross-frontal  mixing  and/or  advection 

13  to  overcome  the  wind  stirring. 

14  (4)  The  thermal  finestructure  along  the  southern  tongue  of  the  Cheju-Yangtze 

15  Front  has  large  variability  locally.  In  the  northeastern  part  of  the  southern  tongue  (south  of 

16  Cheju  Island),  it  is  characterized  by  summer  intensive  thennocline  intrusions  at  the 

17  horizontal  scales  of  tens  kilometers  and  the  vertical  scales  of  a  few  to  40  m.  Similar  but 

18  smaller-scale  haline  finestructure  was  also  observed  in  this  region  (Lee  et  al.,  2003).  These 

19  thermohaline  finestructures  may  be  driven  by  the  vertically-sheared  advection.  In  addition, 

20  Lee  et  al.  (2003)  suggested  that  both  double-diffusive  process  along  isopycnal  surface  and 

21  baroclinity  should  be  taken  into  account.  These  intrusions,  exhibiting  large  Cox  number, 
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1  large  gradient  ratio,  and  large  interquartile-range,  contributes  large  lateral  mixing  related  to 

2  weak  horizontal  density  gradient,  great  cross-frontal  heat  flux,  and  great  lateral  eddy 

3  diffusivity  (Lee  et  al.,  2003;  Park  and  Chu,  2007b). 

4  On  the  other  hand,  the  (inestructure  in/around  the  southwestern  part  of  the  southern 

5  tongue  (southern  Yangtze  Bank)  is  not  as  rich  as  that  in/around  the  northeastern  part  in 

6  warm  seasons.  This  is  attributed  to  weak  cross-frontal  mixing  around  the  front,  which  is 

7  mainly  caused  by  strong  along-frontal  currents,  such  as  the  Taiwan  Warm  Current  and/or 

8  the  Kuroshio  meandering. 

9  (5)  The  profiles  in  the  Tsushima  Warm  Current  branching  region  are  different  even 

10  locally,  and  such  difference  is  more  distinctive  in  the  spring.  The  profile  on  the  branching 

1 1  location  displays  larger  irregular-staircase  thennocline  (at  vertical  scale  of  10-50  m)  in  the 

12  upper  part  (shallower  than  230  m,  like  the  profiles  north  of  the  branching  location)  than  in 

13  the  lower  part  (like  the  profiles  south  of  the  branching  location).  The  strong  horizontal 

14  temperature  gradient  between  the  uplifted  cold  subsurface  Kuroshio  water/the  cold  shelf 

15  water  and  the  warm  Kuroshio  surface/near-surface  water  is  the  source  to  develop  this  large 

16  irregular-staircase.  Related  haline  finestructure,  which  is  more  multilayered,  is  seen  in 

17  several  studies  (Bao  et  al.,  1996;  Lie  et  al.,  1998;  Lee  et  al.,  2003).  The  front  in  this  region 

18  is  partially  density-compensated.  The  lower  bound  of  these  thennohaline  finestructures 

19  (-200  m  depth)  generally  coincides  with  the  depth  of  maximum  salinity  zone  (-34.7  psu) 

20  and  seasonal  pycnocline. 
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Figure  captions 

2  Figure  1.  Various  finestructures:  (a)  a  staircase  (profile  1  observed  at  126. 3°E,  27.0°N),  (b) 

3  a  thermocline  intrusion  (profile  2  at  127. 1°E,  32.8°N;  profile  3  at  126. 3°E,  32.8°N),  (c) 

4  mixed  structures  (profile  4  at  128. 4°E,  34.1°N;  profile  5  at  127.2°E,  33.6°N;  profile  6  at 

5  127. 1°E,  31.8°N).  See  details  in  the  text.  A  thin  black  curve  in  (a)  is  the  vertical 

6  temperature  gradient  of  profile- 1  and  its  scale  is  on  the  top  of  the  figure.  An  inside  figure  of 

7  (a)  is  an  enlargement  of  the  segment  (200-250  m)  of  the  temperature  of  profile- 1. 

8  Figure  2.  (a)  Four  airborne  expendable  bathythermograph  (AXBT)  surveys  conducted  on 

9  18-29  September  1992  (9209),  4-14  February  1993  (9302),  5-14  May  1993  (9305)  and  2- 

10  10  September  1993  (9309).  Contours  indicate  the  bottom  bathymetry  shallower  than  1000 

11  m.  Seven  grey  boxes  indicate  locally-grouped  profiles:  Yellow  Sea  Bottom  Cold  Water 

12  group  (YCWG),  west  of  Cheju  group  (WCG),  south  of  Cheju  group  (SCG),  Cheju- 

13  Tsushima  Front  group  (CTFG),  Yangtze  Bank  group  (YBG),  south  of  Yangtze  Bank  group 

14  (SYBG),  and  Kuroshio  Front  group  (KFG).  (b)  Identified  thermal  fronts  with  a  temperature 

15  distribution  at  25  m  depth  from  the  AXBT  surveys:  Cheju-Yangtze  Front  (CYF),  Cheju- 

16  Tsushima  Front  (CTF),  Tsushima  Front  (TF),  and  Kuroshio  Front  (KF)  (Park  and  Chu, 

17  2007b). 

18  Figure  3.  Temperature  profiles  in  KFG  (location  see  Fig.  2a). 

19  Figure  4.  (a)  A  temperature/salinity  profile  observed  at  127. 14°E  and  28.25°N  in  June  1994, 

20  provided  by  WOD05:  a  back  (gray)  line  for  temperature  (salinity),  (b)  A  vertical 

21  temperature/salinity  gradient  from  the  profiles  of  (a):  a  back  (gray)  line  for  temperature 
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(salinity).  A  horizontal  axis  (a  unit  10"5  m)  was  scaled  by  the  thermal  expansion  and  haline 

2  contraction  coefficients  so  that  the  fluctuations  of  temperature  and  salinity  gradients  of 

3  equal  size  imply  identical  effects  on  density,  (c)  Turner  angle  of  from  the  profiles  of  (a): 

4  see  details  in  the  text.  The  density  ratio  R  corresponding  to  the  Turner  angle  is  marked  at 

5  the  top. 

6  Figure  5.  A  histogram  of  the  Turner  angle  calculated  from  260  profiles  (51-400  m)  in 

7  KFG:  the  profiles  data  were  extracted  from  WOD05  for  a  period  of  1991-2002.  A  bin  size 

8  is  jt/50.  The  density  ratio  R  corresponding  to  the  Turner  angle  is  marked  at  the  top. 

9  Figure  6.  Temperature  profiles  in  CTFG  (see  Fig.  2a). 

10  Figure  7.  Temperature  profiles  in  YCWG  (see  Fig.  2a). 

1 1  Figure  8.  Temperature  profiles  in  WCG  (see  Fig.  2a). 

12  Figure  9.  Temperature  profiles  in  SCG  (see  Fig.  2a). 

13  Figure  10.  Temperature  profiles  in  YBG  (see  Fig.  2a). 

14  Figure  11.  Temperature  profiles  in  SYBG  (see  Fig.  2a). 

15  Figure  12.  Spectra  of  vertical  temperature  gradient  in  the  groups.  KFG  is  divided  into  upper 

16  segment  (upper-KFG,  shallower  than  150m)  and  lower  segment  (lower-KFG,  151-400  m). 

17  Figure  13.  Ensemble-averaged  PDFs  of  vertical  temperature  gradient  in  the  same  groups  as 

18  in  Fig.  10.  Marks  are  placed  at  the  center  of  bin  with  a  size  of  0.05  °C  m'1. 

19  Figure  14.  (a)  A  diagram  of  interquartile-range  vs.  skewness  of  vertical  temperature 

20  gradient  in  the  same  groups  as  in  Fig.  10  (see  details  in  the  text).  The  diagram  is  divided 

21  into  four  zones:  zone-1  (interquartile-range  <  0.05  and  skewness  <  2),  zone-2 
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1  (0.05  <  interquartile-range  <  0.1  and  skewness  <  2),  zone-3  (skewness  >  2),  and  zone-4 

2  (interquartile-range  >0.1  and  skewness  <  2).  (b)  A  histogram  of  profile  population 

3  percentage  of  the  four  zones. 

4  Figure  15.  (a)  A  diagram  of  Cox  number  vs.  gradient  ratio  of  vertical  temperature  gradient 

5  in  the  same  groups  as  in  Fig.  10  (see  details  in  the  text).  The  diagram  has  three  zones:  zone- 

6  a  (Cox  number <0.7  and  gradient  ratio  <  0.1),  zone-)!  (Cox  number  >0.7  and  gradient 

7  ratio  >0.1),  and  zone-y  (0.5  <  Cox  number  <2  and  0.05  <  gradient  ratio  <  0.5).  (b)  A 

8  histogram  of  profile  population  percentage  of  the  three  zones. 

9  Figure  16.  Systematic  of  the  current  system  in  the  East  China  Sea.  The  Kuroshio  following 

10  along  the  continental  shelf  is  separated  into  the  northern  branch,  i.e.  the  Tsushima  Warm 

1 1  Current,  and  the  eastern  northern  branch  at  127-127°E  and  29-30°N,  which  is  assumed  as 

12  the  Tsushima  Warm  Current  branching  location  (Huh,  1982;  Lie  et  al.,  1998;  Hsueh,  2000; 

13  Furey  and  Bower,  2005).  The  Taiwan  Warm  Current,  which  passes  through  the  Taiwan 

14  Strait  and  then  flows  over  the  continental  shelf,  contributes  to  the  Tsushima  Warm  Current 

15  in  warm  seasons  (Isobe,  1999;  Ichikawa  and  Beardsley,  2002).  The  AXBT  stations  are 

16  marked  by  a  cross.  Three  individual  AXBT  profiles  marked  by  an  open  circle  are  used  for 

17  calculating  spectral  and  statistical  values  in  the  Tsushima  Warm  Current  branching  region: 

18  a  station  at  the  branching  location  labeled  by  “A”  and  the  two  neighboring  stations  labeled 

19  by  “B”  and  “C”  (Stn-A,  Stn-B,  and  Stn-C,  respectively).  A  box  is  for  Fig.  17. 

20  Figure  17.  Temperature  profiles  in  the  Tsushima  Warm  Current  branching  region  (a  box  in 

21  Fig.  16). 
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Figure  18.  Temperature  profiles  (a),  vertical  temperature  gradient  (b),  cumulative  variance 

2  of  vertical  temperature  gradient  (c),  and  distribution  probability  of  vertical  temperature 

3  gradient  (d)  of  three  stations  in  the  Tsushima  Wann  Current  branching  region  (see  Fig.  16; 

4  Stn-A  for  the  first  row,  Stn-B  for  second,  and  Stn-C  for  third).  In  (b)  and  (c)  9209  is  correct 

5  for  scale  at  top  and  remaining  curves  are  offset  by  0.5  °C  m'  and  1  °C  m'  ,  respectively. 

6  Figure  19.  Spectra  of  vertical  temperature  gradient  of  the  three  stations  in  the  Tsushima 

7  Warm  Current  branching  region  (see  Fig.  16). 
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1  Table  1 .  Cox  number  averaged  over  the  group. 
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Figure  1.  Various  finestructures:  (a)  a  staircase  (profile  1  observed  at  126. 3°E,  27.0°N),  (b)  a  thermocline 
intrusion  (profile  2  at  127. 1°E,  32.8°N;  profile  3  at  126. 3°E,  32.8°N),  (c)  mixed  structures  (profile  4  at 
128.4°E,  34.1°N;  profile  5  at  127.2°E,  33.6°N;  profile  6  at  127.1°E,  31.8°N).  See  details  in  the  text.  A 
thin  black  curve  in  (a)  is  the  vertical  temperature  gradient  of  profile  1  and  its  scale  is  on  the  top  of  the 
figure.  An  inside  figure  of  (a)  is  an  enlargement  of  the  segment  (200-250  m)  of  the  temperature  of 
profile  1. 
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Figure  2.  (a)  Four  airborne  expendable  bathythermograph  (AXBT)  surveys  conducted  on  18-29  September 
1992  (9209),  4-14  February  1993  (9302),  5-14  May  1993  (9305)  and  2-10  September  1993  (9309). 
Contours  indicate  the  bottom  bathymetry  shallower  than  1000  m.  Seven  grey  boxes  indicate  locally- 
grouped  profdes:  Yellow  Sea  Bottom  Cold  Water  group  (YCWG),  west  of  Cheju  group  (WCG),  south  of 
Cheju  group  (SCG),  Cheju-Tsushima  Front  group  (CTFG),  Yangtze  Bank  group  (YBG),  south  of 
Yangtze  Bank  group  (SYBG),  and  Kuroshio  Front  group  (KFG).  (b)  Identified  thermal  fronts  with  a 
temperature  distribution  at  25  m  depth  from  the  AXBT  surveys:  Cheju-Yangtze  Front  (CYF),  Cheju- 
Tsushima  Front  (CTF),  Tsushima  Front  (TF),  and  Kuroshio  Front  (KF)  (Park  and  Chu,  2007b). 


42 


1 

2 

3 

4 

5 

6 


Figure  4.  (a)  Temperature  (black  curve)  and  salinity  (grey  curve)  profiles  observed  at 
127.14°E  and  28.25°N  in  June  1994  (from  WOD05),  (b)  vertical  temperature  gradient 
multiplied  by  the  thermal  expansion  coefficient  ( a )  (black  curve)  and  haline  gradient 
multiplied  by  the  salinity  contraction  coefficient  (grey  curve)  (unit:  10'5  irf '),  and  (c) 
Turner  angle  and  corresponding  density  ratio  R  (unit  marked  at  the  top). 
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2  Figure  5.  Flistogram  of  the  Turner  angle  (or  density  ratio  R)  calculated  from  260  profiles 

3  (5 1-400  m)  in  the  KFG  with  a  bin  size  of  jt/50.  The  data  are  extracted  from  WOD05  for  the 

4  period  of  1991-2002.  The  density  ratio  R  is  marked  at  the  top. 
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3  Figure  10.  Temperature  profiles  in  YBG  (location  see  Fig.  2a). 
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Figure  11.  Temperature  profiles  in  SYBG  (location  see  Fig.  2a). 
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3  Figure  12.  Spectra  of  vertical  temperature  gradient  in  various  groups  with  KFG  divided  into 

4  upper  segment  (upper-KFG,  shallower  than  150  m)  and  lower  segment  (lower-KFG,  151  — 

5  400  m). 
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Figure  13.  Ensemble-averaged  PDFs  of  vertical  temperature  gradient  with  bin  size  of  0.05 
°C  in1  for:  (a)  upper-KFG,  (b)  lower-KFG,  (c)  CTFG,  and  (d)  SCG. 
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2  Figure  14.  (a)  Scatter  diagram  of  interquartile -range  vs.  skewness  of  vertical  temperature 

3  gradient  for  the  same  groups  in  Fig.  10  (see  details  in  the  text)  with  four  zones:  zone-1 

4  (interquartile-range  <0.05  and  skewness  <  2),  zone-2  (0.05  <  interquartile-range  <0.1 

5  and  skewness  <  2),  zone-3  (skewness  >2),  and  zone -4  (interquartile-range  >0.1  and 

6  skewness  <  2).  (b)  Histograms  of  occurrence  percentage  of  the  four  zones. 
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2  Figure  15.  (a)  Scatter  diagram  of  Cox  number  vs.  gradient  ratio  of  vertical  temperature 

3  gradient  for  the  same  groups  as  in  Fig.  10  (see  details  in  the  text)  with  three  zones: 

4  zone-a  (Cox  number<0.7  and  gradient  ratio<0.1),  zone-P  (Cox  number>0.7  and 

5  gradient  ratio>0.1),  and  zone-y  (0.5<Cox  number<2  and  0.05 <gradient  ratio<0.5). 

6  (b)  Histograms  of  occurrence  percentage  of  the  three  zones. 
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Figure  16.  Systematic  of  the  current  system  in  the  East  China  Sea.  The  Kuroshio  following 
along  the  continental  shelf  is  separated  into  the  northern  branch,  i.e.  the  Tsushima 
Warm  Current,  and  the  eastern  northern  branch  at  127-127°E  and  29-30°N,  which  is 
assumed  as  the  Tsushima  Warm  Current  branching  location  (Huh,  1982;  Lie  et  ah, 
1998;  Hsueh,  2000;  Furey  and  Bower,  2005).  The  Taiwan  Warm  Current,  which 
passes  through  the  Taiwan  Strait  and  then  flows  over  the  continental  shelf,  contributes 
to  the  Tsushima  Warm  Current  in  warm  seasons  (Isobe,  1999;  Ichikawa  and  Beardsley, 
2002).  The  AXBT  stations  are  marked  by  a  cross.  Three  individual  AXBT  profiles 
marked  by  an  open  circle  are  used  for  calculating  spectral  and  statistical  values  in  the 
Tsushima  Warm  Current  branching  region:  a  station  at  the  branching  location  labeled 
by  “A”  and  the  two  neighboring  stations  labeled  by  “B”  and  “C”  (Stn-A,  Stn-B,  and 
Stn-C,  respectively).  A  box  is  for  Fig.  17. 
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3  Figure  18.  (a)  Temperature  profiles,  (b)  vertical  temperature  gradient,  (c)  cumulative 

4  variance  of  vertical  temperature  gradient,  and  (d)  probability  distribution  function  of 

5  vertical  temperature  gradient  for  three  stations  in  the  Tsushima  Warm  Current 

6  branching  region  (location  see  Fig.  16)  with  the  first  row  for  Stn-A,  second  row  for 

7  Stn-B,  and  third  row  for  Stn-C.  In  (b)  and  (c)  the  profile  9209  is  in  correct  scale  shown 

8  at  the  top  and  remaining  profiles  are  offset  by  0.5  °C  m'  and  1  °C“  m'  ,  respectively. 
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3  Figure  19.  Spectra  of  vertical  temperature  gradient  of  the  three  stations  in  the  Tsushima 

4  Warm  Current  branching  region  (location  see  Fig.  16). 
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